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INTRODUCTION 


For  some  time  there  has  been  considerable  concern  with  average  cost 

computations  in  inventory  systems  with  non-unit  order  sizes  and  variable 

lead  times  when  following  a stationary  (S,s)  policy.  In  this  paper  the 

derivation  of  the  stationary  distribution  of  the  inventory  position  was 

motivated  to  a large  extent  by  Izzet  Sahin's  work  reported  in  [2];  however, 

» 

in  that  paper  much  time  was  spent  on  deriving  finite  horizon  equations 
and  then  taking  limits.  We  circumvent  the  latter  by  applying  the  well 
known  age  distribution  for  renewal  processes.  The  work  also  generalizes 
Silvers  [3]  results  on  the  stationary  distribution  of  inventory  position. 

In  the  first  chapter  we  derive  the  stationary  distribution  of  the 
inventory  position  and,  under  the  assumption  of  at  most  one  order  out- 
standing, derive  the  joint  stationary  distribution  of  the  inventory 
position  and  the  number  of  orders  outstanding.  The  latter  distribution 
pernits  us  to  compute  approximately  the  average  cost  for  given  demand  and 
lead  time  distributions. 

In  the  second  chapter  we  analyze  an  example  with  a specific  form  of 
demand  and  lead  time  distribution;  developing  recursion  formulae  for 
fast  computation  of  the  joint  stationary  distribution  of  inventory  position 
and  number  of  orders  outstanding.  Also  given  is  a listing  of  a computer 
program  which  implements  these  recursion  formulae  and  finds  the  minimum 
cost  value  of  s for  a fixed  S. 

The  third  chapter  takes  a heuristic  approach  and  gives  an  average 
cost  equation  based  on  this  author's  intuitive  notions  on  the  ergodic 
properties  of  some  of  the  processes  involved,  plus  some  outright  approxima- 
tions. No  claims  to  rigor  are  made  here,  but  then  of  course  the  analysis 
is  short  and  sweet  and  its  adequacy  can  be  tested  via  simulation. 
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CHAPTER  I 


RIGOROUS  ANALYSIS 

Consider  a continuous  review  inventory  system.  Requisitions  for 
inventory  items  occur  at  random  points  in  time  and  the  size  of  each 
requisition  is  also  random.  The  ordering  policy  is  (S,s)  and  the  lead 
times  are  random. 

If  I (t)  is  the  inventory  position  (on  hand  plus  on  order)  at  time  t 
then  the  times  at  which  orders  are  placed  are  renewal  points  for  the 
Ip(t)  process.  Define 

N(t)  = number  of  orders  placed  by  time  t 

T = time  between  the  n*"'1  and  n-lSt  order 
n 

N(t) 

and  V(t)  = t - £ T 

i n 

n=l 

V(t)  is  the  time  since  the  last  order  was  placed  and  is  sometimes  called 
the  age  at  time  t. 

If  we  let  an  asterisk  denote  the  stationary  version  of  a process 
then  it  is  well  known  that  (see  [3]  pp.  45) 

Fv*(y>  = ET*  ° [i"FT*(u)]du 

so  if  A = S-s  then  for  s < x <_  S and  D(*)  cumulative  demand 
P[I*  = X]  = f P[I*  = X | V*  ->  y]  fy*  (y)dy 

= ^ /“  P[D(y)  « S - X|  V*  = y]  [1-F  *(y)]  dy 
ET  T 


f P[D(y) 


S - X | D(y)  < A]  [1-F  *(y)]dv 

T 


f Pln(y) 


S - X] dy 
3 


(1) 


p 


since  V =y<  = > D(y)  < A,  and  a renewal  at  -y 
Moreover  /°°  P[D(y)  - S - X]dy 

= E fo  I[D(y)  - S - X]dy 


(2) 


» E [amount  of  time  cumulative  demand  was  at  S-X] 

P [cumulative  demand  was  ever  at  S-X] 

X 

where  X 1/E  (time  between  requisitions). 

Let 

B(X)  = CDF  of  the  requisition  size  and  B^(X)  its  n-fcld 

convolution  (B  (X)  » 1) 
o 

A(X)  = CDF  of  cumulative  requisitions  and  A (X)  its  n-fold 

n 

convolution  (A  (X)  » 1) 
o 

CO 

H (X)  = T.  B (X)  ; H (o)  = o 

i n 

n=l 

It  is  well  known  that  H is  the  renewal  function  (expected  cumulative 
renewals)  for  the  process  with  inter-renewal  distribution  B,  i.e.,  the 
time-between-renewals  distribution  is  just  the  distribution  of  requisition 
size  and  H(X)  is  the  expected  number  of  requisitions  in  a total  demand  of  X. 


h(X)  H H (X)  - H(X-l)  = P [Renewal  at  X]  x > o 


and  a renewal  at  X for  this  process  is  the  same  as  the  event  that  cumulative 
demand  is  ever  at  X,  hence  from  (2) 


h(X) 


X / P[D(y) 


X]  dy 


(3) 


; 


4 


Since  the  T are  lid  T and  T have  the  sane  distribution  and  it  is  clear 
n n 

that 

T > y <=>  U(y)  < A 


Hence 


F.T  = / I»[T  > v]dy  - / P[D(y)  < A]  dy 

o ‘ o 


1 


A-l 

OO  00 

= l I P[D(y>*=  X]  dy  + / P[D(y)  = o]  dy 
X-l  ° 

A-l 

= V 

x-l 


h(X) 

X 


+ / P[D(v)  = o]  dy 
o 


= 1/X  [H(A-l)  - H(o) ] + 1/X 


= y [H(A-l)  + 1] 


(4) 


Use  the  natural  extension  of  (3)  to  get  h(o) 
equations  1-4  to  get 


X / P[D(y) 
o 


ol  dy  = 1 and 


P l I = X]  = 


. h(S-X) 


1 + H(A-l) 


s < X < S 


Let 


L = lead  time  of  an  order 
Z(t)  = number  of  orders  outstanding  at  time  t 


Suppose  the  probability  that  the  lead  time  demand  is  greater  than  A-l 
is  negligible.  Then  we  can  safely  assume  that  there  is  never  more  than  one 
order  outstanding,  i.e.,  Z(t)  is  o or  1 almost  surely.  It  then  follows 
that 

P[I*  = X,  Z*  = o]  - f P[D(y)  = S - X]  P [ L y]  dy 


I 


1+H(A-1) 


/ P [ D ( v ) = S - X]  P [L  < y [dy 
o — 


and 


* * 

P[ID  = X,  z 


= 1] 


h(S-X) 

1+H(A-1, 

5 


P[I*  = X,  Z*  = o] 


since  P [ I * = X]  = P[I*  « X,  Z*  - 1]  + P 1 1*  - X,  Z*  = o] 


these  last  equations  can  be  used  to  compute  approximately 
i:  [Average  backorder  anu  holding  cost] 


. * * + * * - 
= HC  • E(lp  - Z A)  + HB  • E (Ip  - Z A) 

*■  ■ 

where  we  have  approximated  the  outstanding  order  size  by  A (it  is  in  general 
larger)  and 

HC  = holding  cost  per  item  per  unit  time 
HB  = backorder  cost  per  item  per  unit  time 
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CHAPTER  II 


EXAMPLE 

If  the  time  between  requisitions  is  exponential  with  parameter  A 
ami  the  requisition  sizes  are  geometric  with  parameter  p,  i.e.. 


p(x)  = pq 


x * 1,2,...  ; q = 1 - p 


P[D(y)  = x]  = 


( ~>.y 

e 


-Ay  x 

e q 


v i_ 

j-iV-V  j! 


B (x)  = 
n 


v ( k-1  \ n k-n 

u W p q 


X 

x x /k-1 

l 

njx) 

= Z Z I 

, \n-l 

n=l  k^n  ' 

n=l 

k / 

k_1\  n k-n 

i n-1  ) P q 

= 1 

n=l  ' 

n k-n 
P q 


* k 1 /k-l\  n k-1  -n 

= s p - ) p q 

k=l  n=o  \ n / 


p = px 


ET  = 1/A  (1  + p(A-l)) 
( 1 x - 


o 


7 


k 


and 


‘ P[D(y)  = x Jdy 


-Az 


e ’ / A)  Z a (x, z)  x > o 


J-l 


J 


wnere 


a 


aj<x’z> 


k, 


(p/q> - v (Az)  kl 


K=0 


Suppose  the  support  of  L is  a finite  lattice,  i.e.,  L assumes  only 
a finite  number  of  values.  Let 


P[L  = gj  = y± 


1 - 1,2, ...m 


Then  we  can  write 


* A 

P [ I = x,  7.  = o' 

P 


/ P[D(y)  = x]  P[L  < y]dy 


St 

y^^  / [D Cy ) *=  x]dy  + y -ry0)  / ' P[D(y)=x]dy 


+ ...  + (y1  + ...  + ym)  r P[D(y)  = x]dy 


= yx  f P [D (y)  * x]dy  + yn  /“  P[D(y)  = x]dy 


+ ...  + y^  / P [D(y)  = x]dy 


8m 


ra 


t ye  gi/A 
i=l  1 


x = o 


m x 

(qXe  8i/X)  z a (x,g  ) 
i-1  1=1  J 1 


x > o 
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For  computations  let 


J k 

b •=  T.  ( A 2 ) / k ! 

J k=o 

and 

c ^ = (Az) ^ / j ! 
so  that 

= a.  [1  + Azc^.  / (j+1)  bj  ] [“^-]  (p/q) 


This  recursion  affords  a straight  forward  and  fast  algorithm  for  computing 
* * 

P[I  = x,  Z = o]  . 

P 


Remark:  For  this  example  empirical  investigations  show  that  the  average  cost 

is  a convex  function  of  s. 

Appendix  A contains  a listing  of  a computer  program  which  implements 
these  recursion  formulae  and  finds  the  minimum  cost  value  of  s for  a given 
order-up-to  point,  backorder  ost,  holding  cost,  requisition  rate,  and 
average  requisition  size.  To  use  it 


CALL  OPTRPT  (S,IS,  BC,  IIC,  CO,  L,  P) 


INPUT:  BC  = unit  backorder  cost 

HC  = unit  holding  cost 
S = order  up  to  point 
L requisition  rate 

P = recipricol  of  the  average  requisition  size 

OUTPUT:  IS  = optimal  reorder  point 

CO  = optimal  cost 


i 


1 

; i 

j i 


I 
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CHAPTER  III 
HEURISTIC  ANALYSIS 

As  an  alternative  approach  consider  points  in  tine  at  which  orders 
arrive.  The  average  cost  incurred  between  two  of  these  points  is  given 
by  (see  figure  below) 

M m 

CM(M)  = CH  M+l  Z (N  " nEr)+  + CB  " (N  - nEr)“ 

n=o  n=o 

Lo 

M M 

= [CH  I (N-nEr)  + CB  I (N-nEr) “] 
n=o  n=o 

wiiere  Lo  = time  between  the  arrival  of  the  two  orders 
M = number  of  requisitions  in  Lo 
Ei  = expected  requisition  size 
N * net  inventory  just  after  first  order  arrives 
Cg  = unit  backorder  cost 
* unit  holding  cost 
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1 

I 

i 


Let  ET  be  the  expected  time  between  the  placement  of  orders.  Since 
orders  are  being  placed  at  a rate  1/ET,  so  must  the  orders  be  arriving  at 
the  same  rate. 

So  lets  be  bold  and  replace  ELo  with  ET  to  get  EM  - AELo  - AET . It 
looks  like  the  necessary  ergodic  properties  should  hold  so  that  we  can 
write 


ecm(n) 


IChE 


M 
E 

n“o 


(N  - nEr)  + C 


B 


M 

E Z 


(N  - nEr)"] 


n“o 


AET+1 


Assuming  there  is  always  more  than  one  order  outstanding  the  order  of 
smoothing  in  the  numerator  is  arbitrt ry  im  e and  N are  independent. 
However,  M depends  on  Lo  inspite  of  the  notation.  Replace  Lo  with  T to  get 
the  distribution  of  M. 

For  the  distribution  of  N (net  inventory  at  the  time  that  an  order 
arrives)  note  that  if  orders  cannot  pass  then 


.V  - S - D(L) 


Because  in  this  case  all  orders  outstanding  at  t-L  are  in  at  time  t,  all 
orders  placed  after  time  t-L  are  still  out  at  time  t and  the  inventory 
position  at  t-L  is  S since  t-L  is  an  order  point  if  t is  the  time  the  order 
arrived.  Hence 

P [N  - x]  - P [S  - D(L)  = x] 

= P[D(L)  - S - x] 

If  L is  on  a finite  lattice  then  the  above  is  not  difficult  to  compute 
if  we  have  the  distribution  of  D(y);  for  example,  the  cumulative  demand 
distribution  derived  in  the  analytic  approach  of  the  first  chapter. 


1 


12 


REFERENCES 


E.A.  Silver,  Chung-Mei  ilo,  Robert  Deemer,  "Cost-Minimizing  Inventory 
Control  of  Items  Having  a Special  Type  of  Erratic  Demand  Pattern," 
Infor.,  Vol  9,  No.  3,  November  1971,  198-219. 

I.  Sahin,  "On  the  Stationary  Analysis  of  (S,s)  Inventory  Systems 
With  Lead  Times,"  Working  Paper  76-10,  Faculty  of  Management  Sciences, 
University  of  Ottawa. 

S.  Ross,  Applied  Probability  Models  with  Optimization  Applications, 
Holden-Day,  San  Francisco,  1970. 


DISTRIBUTION 


I 


i 


i 


i 

i 


COPIES 

1 Deputy  Under  Sec'y  of  the  Army,  ATTN:  Office  of  Op  Resch 

Headquarters,  US  Army  Materiel  Development  6 Readiness  Command 
1 DRCPA-S 

1 DRCMM-RS 

2 Defense  Logistics  Studies  Info  Exchange,  DRXMC-D 

2 Commander,  US  Army  Logi^ics  Center,  Ft.  Lee,  VA  23801 

1 Commander,  US  Army  Tank-Automotive  Research  & Development  Cmd, 

ATTN:  DRDTA-V , Warren,  MI  48090 
1 Commander,  US  Army  Tank-Automotive  Materiel  Readiness  Cmd, 

ATTN:  DRSTA-S , Warren,  MI  48090 
1 Commander,  USA  Armament  Command,  Rock  Island,  ILL  61201 

1 Commander,  USA  Electronics  Command,  Ft.  Monmouth,  NJ  07703 

1 Commander,  USA  Missile  Research  4 Development  Cmd,  Redstone 

Arsenal,  Ala  3580) 

1 Commander,  USA  Missile  Readiness  Command,  Redstone,  Ars,  Ala  35809 

1 Commander,  USA  Troop  Support  Cmd,  4300  Goodfellow  Blvd. , St. 

Louis,  MO  63120 

1 Commander,  USA  Aviation  Systems  Cmd,  P.O.  Box  209,  St.  Louis, 

MO  63120 

1 Commander,  Army  Automated  Logis.ic:  Mgt  Systems  Agency,  Box 

14505,  St.  Louis,  '10  63168 

1 Director,  DARCOM  Logistics  Systems  Support  Agency,  Letterkenny 

Army  Depot,  Chambersburg , PA  17201 
1 Commander,  Maintenance  Mgt  Center,  Lexington-Blue  Grass  Army 

Depot,  Lexington,  KY  40507 

1 Director,  Army  Management  Engineering  Training  Agency,  Rock 

Island  Arsenal,  Rock  Island,  ILL  61202 
1 Commandant,  US  Army  Logistics  Mgt  Center,  Ft.  Lee,  VA  23801 

1 Defense  Logistics  Agency,  Cameron  Sta,  Alexandria,  VA  22314 

4 Dep  Chf  of  Staff  (I&L),  HQ  USMC-LMP-2,  ATTN:  MAJ  Sonneborn,  Jr., 

Wash.,  DC  20380 

10  Defense  Documentation  Center,  Cameron  Sta.,  Alexandria,  VA 

1 Commander,  US  Air  Force  Logistics  Command,  WPAFB,  Dayton, 

Ohio,  ATTN:  AFLC/XRS  45433 

1 US  Navy  Fleet  Materiel  Support  Office,  Naval  Support  Depot, 

Mechanicsburg , PA  17055 

1 Mr.  James  Prichard,  Navy  Supply  Systems  Cmd,  Dept  of  US  Navy, 

Wash. , DC 

1 George  Washington  University,  Inst  of  Management  Science  & 

Engineering,  707  22nd  St.,  N.W.,  Washington,  DC 
1 Naval  Postgraduate  School,  ATTN:  Dept  of  Opns  Anal,  Monterey, 

Calif  93940 

1 Air  Force  Institute  of  Technology,  ATTN:  SLGQ,  Head  Quantita- 

tive Studies  Dept.,  Dayton,  Ohio  43433 


' 


16 


The  Army  Library,  Room  1A518,  Pentagon,  Wash.,  DC  20310 
US  Army  Military  Academy,  West  Point,  NY 

Logistics  Management  Institute,  4701  Sangamore  Road,  Wash. , DC  200i6 
University  of  Florida,  ATTN:  Dept  of  Industrial  & Systems 

Engineering,  Gainesville,  Fla 
RAND  Corp.,  1700  Main  St.,  Santa  Monica,  CA  90406 
Office,  Asst  Sec 1 y of  Defense  (Inst  6 Logistics)  ATTN:  Mr.  George 
Minter,  Pentagon,  Wash.,  DC  20310 
US  Army  Materiel  Systems  Analysis  Activity,  ATTN:  AMXSY-CL, 

Aberdeen  Proving  Ground,  MD  21003 
Commander,  US  Army  Logistics  Center,  ATTN:  Studies  Analysis 

Div.,  Concepts  & Doctrine  Directorate,  Ft.  Lee,  VA  23801 
ALOG  Magazine,  ATTN:  Tom  Johnson,  USALMC,  Ft.  Lee,  VA  23801 

Commander,  Air  Force  Logistics  Crad,  ATTN:  AFLC/AQMLE,  WPAFB , 

Dayton,  Ohio  43433 

Operations  Research  Center,  3115  Etcheverry  Hall,  University 
of  California,  Berkeley,  CA  94720 
Dr.  Jack  Muckstadt,  Dept  of  Industrial  Engineering  A Operations 
Research,  Upson  Hall,  Cornell  University,  Ithaca,  NY  14890 
Prof  Herbert  P.  Oalliher,  Dept  of  Industrial  Engineering, 

University  of  Michigan,  Ann  Arbor,  MI  48104 
Mr.  Ellwood  Hurford,  Scientific  Advisor,  Army  Logistics  Center, 

Ft.  Lee,  VA  23801 

Prof  Robert  '!.  Stark,  Dept  of  Stat  A Computer  Sciences, 

University  of  Delaware,  Newark,  DEL  19711 
Prof  E.  Gerald  Hurst,  Jr.,  Dept  of  Decision  Science,  The  Wharton 
School,  University  of  Penna.,  Phila.,  PA  19174 
Logistics  Studies  Office,  DRXMC-LSO,  ALMC , Ft.  Lee,  VA  23801 
Procurement  Research  Office,  DRXMC-PRO,  ALMC,  Ft.  Lee,  VA  23801 
Commander,  ' ’ Army  Aviation  Systems  Cmd,  ATTN:  DRSAV-D,  P.O. 

Box  209,  St.  Louis,  MO  63166 

Commander,  US  Army  Electronics  Cmd,  ATTN:  DRSEL-SA,  Ft.  Monmouth, 

NJ  07703 

Commander,  US  Army  Test  A Evaluation  Cr,.d,  ATTN:  DRSTE-SY’, 

Aberdeen  Proving  Ground,  MD  21005 
Commander,  US  Army  Armament  Cmd,  ATTN:  DRSAR-SA,  Rock  Island, 

IL  61201 

Commander,  US  Army  Troop  Support  Command,  ATTN:  DRSTS-G, 

4300  Goodfellow  Blvd,  St.  Louis,  MO  63120 
Commander,  US  Army  Mobility  Equipment  Research  A Development 
Cmd,  ATTN:  DRXFB-0 , Ft.  Belvoir,  VA  22060 

Commander,  US  Army  Natick  Research  A Development  Cmd,  ATTN: 

DRXNM-0,  Natick,  MA  01760 

Commander,  USA  Missile  Cmd,  ATTN:  DRSMI-C,  Redstone  Ars,  AL  35809 
Commander,  USA  Missile  Cmd,  ATTN:  DRSMI-D,  Redstone  ARS,  AL  35809 
Commander,  USA  Armaments  Research  A Development  Cmd,  Implementation 
Task  Force,  ATTN:  DRCSA-RP-SE  (Mr.  Larry  Ostuni)  Picatinny 
Arsenal,  Dover,  NJ  07801 


17 


